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ABSTRACT 

We use measurements from the Planck satellite mission and galaxy redshift surveys over 
the last decade to test three of the basic assumptions of the standard model of cosmology, 
ACDM: the spatial curvature of the universe, the nature of dark energy and the laws of gravity 
on large scales. We obtain improved constraints on several scenarios that violate one or more 
of these assumptions. We measure wo = —0.94 ±0.17 (18% measurement) and 1 ± rca = 
1.16±0.36(31% measurement) for models with a time-dependen t equation of state, which is 
an improvement over current best constraints (lAubourg et al.l201^ . In the context of modified 
gravity, we consider popular scalar tensor models as well as a parametrization of the growth 
factor. In the case of one-parameter f{R) gravity models with a ACDM background, we 
constrain < 1>36 x 10“^ (la C.L.), which is an improvement by a factor of 4 on the 
current best (IXull2Q15k . We provide the very first constraint on the coupling parameters of 
general scalar-tensor theory and stringent constraint on the only free coupling parameter of 
Chameleon models. We also derive constraints on extended Cham eleon models, impr oving 
the constraint on the coupling by a factor of 6 on the current best (iHojiati et al.ll2Ql ll) . We 
also measure 7 = 0.612 ± 0.072 (11.7% measurement) for growth index parametrization. We 
improve all the current constraints by combining results from various galaxy redshift surveys 
in a coherent way, which includes a careful treatment of scale-dependence introduced by 
modified gravity. 

Key words: gravitation; modified gravity; galaxies: statistics; cosmological parameters; 
large-scale structure 


1 INTRODUCTION 


Since its development a century ago, General Relativity (GR) has 
eonsistently provided a ver y successful framework to describe th e 
evolution of our Universe |Peeblesiri98Ql : IPavis & Peebles! 1 1983h . 
Nowadays, the prediction of GR for the growth of the large scale 
strueture that we observe around us, is reaching great precision 
as cosmic microwave background (CMB) measurements are pro¬ 
viding us with impressively accurate estimates o f the cosmologi- 
eal parameters (jPlanck Collaboration et al.ll2015ah . Yet, the excite¬ 
ment about the advances of observational cosmology is accompa¬ 
nied by the awareness that we face some major challenges. While 
the standard cosmological model, based on the laws of GR, pro¬ 
vides a very good fit to existing data, it relies on a universe of 
which we understand only ~5% of the content. The remaining en¬ 
ergy budget comes in the form of dark matter (~27%), respon¬ 
sible fo r the clustering of structure, and the cosmological con¬ 
stant A (lEinsteinlll915h (~ 68%), responsible for the phase of ac¬ 
celerated expansion recently entered by the universe. In particu¬ 
lar, the physical understanding of cosmic acceleration represents 


one of the most important challenges in front of modern physics. 
While A is in good agreeme nt with available data, e .g. baryon 
acous ti c osci ll ations (BAG) ([Eisenstein et al. ItOO^: Cole et al.l 
200^1 iHutsil I 2 OO 6 I: lKazinet al.1 I201QI: IPercival et all l201Ql : 


Reid et al.l |201Q|: lAubourg et al.l l20l4 lAnderson et al.l l2014a|3^ . 


Goobar & Leibundsut 1201 ll; ISuzuki et al.l 20121; Rodnev et al.l 

2014h. and CMB ( 

Planck Collaboration et al.l d2015ab. WMAP9 

Bennett et al.l d2013 

)) observations, it suffers from the coincidence 


and fine tuning problems dWeinberd Il989l : ICarro iil l200lh . Sev- 
eral alternatives to A have been proposed in the two decades 
since the discovery o f cosmic acceleration dRiess et al.l 1 19981 : 
IPerlmutter et al.l Il999h . and they can be roughly divided into 
two classes. The first class, to which we will refer as modified 
gravity (MG), corresponds to modifications of the laws of grav¬ 
ity on large scales, designed to achi eve self accelerating solu¬ 
tions when matter becomes negligible dSilvestri & Troddenll2009l : 
IClifton et~^ l2012h : alternatively, one can introduce a dynami¬ 
cal degre e of freedom, commonly dubbed dark energy (DE; first 
coined bv iHuterer & Turned dl999lfi . which is smoothly distributed 
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and starts to dominate the evolution of the Universe at late times 
(ICoDeland et al.ll2006h . 

Undoubtedly, one of the important tasks for modern cosmol- 
ogists, is to perform precision tests of the standard model of cos¬ 
mology (ACDM) and identify areas of tension. In a joint effort, 
one needs also to explore the parameter space of alternative mod¬ 
els. Even though with the current constraints from data, any depar¬ 
ture from ACDM is likely to be small and challenging to detect, 
we are in a unique position to test GR, and the other assumptions 
of ACDM, to unprecedented precision with modern observational 
probes. The three basic assumptions of ACDM which are popularly 
tested are the curvature of the universe, the nature of dark energy 
and the laws of gravitational interaction on large scales. The curva¬ 
ture of the universe can be explored by allowing a curvature density 
parameter, Qk, to be different from zero and free to vary. As for 
the nature of dark energy, we will focus on smoothly distributed 
models where it suffices to test for the deviation of the equation 
of state parameter, w, from -1, which is the value it assumes if 
the acceleration is driven by A. We will consider both a constant 
w as well as a time-dependent one, resorting to the popular CPU 

i )arametrization in terms of wq and u;a, i.e. w — wq -\- Waj^ 
Chevallier & Polarskil200ll : lLindeil2003h . Finally, we will explore 
the nature of gravity by replacing GR with various modified gravity 
models, including Chameleon-type scalar-tensor theories and pop¬ 
ular parametrizations of the growth rate. All these alternatives that 
we consider in our analysis, affect, in one way or another, the rate at 
which large scale structures grow. Models of smoothly distributed 
dark energy, which does not cluster, modify only the background 
dynamics of the universe, but this still has an impact on the rate 
at which structure forms. On the other hand, models of modified 
gravity generally modify both the background and perturbation dy¬ 
namics, leading to a significant effect on the growth rate. 

Modern galaxy redshift surveys, have successfully measured 
the growth ra te using Redshift Space Distortions (hereafter RSD ; 
iKaisej (Il987h ). which is the distortion induced in the galaxy corre¬ 
lation function by the peculiar velocity component of the galaxy 
redshift. Hence, on linear scales, RSD offers a handle both on 
the distribution of matter over-density and peculiar velocity of 
galaxies. Recent galaxy redshift surveys have provided the mea¬ 
surement of fas{z) up to redshift z = 0.8 , where / is the 
growth rate, i.e. the logarithmic derivative of the growth factor, 
and cTs is the rms amplitude of matter fluctuations in a sphere 
of radius 8 h“^Mpc. In this paper, we will test all the three as- 



veys as they have different assumptions. We have looked into these 
assumptions and possible systematic while combining results from 
the different survey and also proposed a way to test scale depen¬ 
dence for modified gravity models using these results. 


2 THEORY 

In exploring the power of RSD data to constrain deviations from 
the standard cosmological scenario, we consider several alterna¬ 
tive models, divided into dark energy models that modify the back¬ 
ground expansion history without introducing any clustering de¬ 


gree of freedom, and those that instead modify only the dynam¬ 
ics of perturbations while keeping the background fixed to ACDM. 
In the former case we consider one and two parameter extensions 
of the standard scenario, corresponding to different equations of 
state for dark energy or a non zero spatial curvature. More specif¬ 
ically we consider: a u;CDM universe, where the equation of state 
for dark energy is a constant parameter that can differ from the 
ACDM value w — —1; a (u;o,rr’a)CDM universe, in which the 
equation of state for dark energy is a function of time and is 
approximation to exact solutions of the scalar field equation of 
motion, i.e. the Chevallier-Polarski-Linder (CPU) parameterization 
w = u;o + u;a(l — a); a oACDM universe which can have a spatial 
curvature different from zero, parameterized in terms of the corre¬ 
sponding fractional energy density Q.k- In the case of models that 
modify the equations for the evolution of perturbations, we analyze 
Chameleon-type scalar-tensor theories, f{R) gravity and a time de¬ 
pendent parametrization of the growth rate. 

We use the publicly av ailable Einstein-Boltzmann solver MG- 
CAMB llHoiiatietal .11201 llfl to evolve the dynamics of scalar per¬ 
turbations and obtain predictions to fit to our data set for all the 
models considered, except for the rr’a)CDM cas e. This lat¬ 
ter ne eds to be treated instead through the PPF module (iFang et al.l 
l2008l) in CAMB 0. While the implementation of the non-clustering 
dark energy models is trivial, in the following we shall describe in 
more detail the implementation of the modified gravity models. 


2.1 Scalar-tensor theories 


Going beyond simple extensions of the standard model and non¬ 
clustering dark energy models, one needs to take into consideration 
also the modifications to the equations for cosmological perturba¬ 
tions. Given the cosmological probes that we consider in our anal¬ 
ysis, it suffices for us to focus on linear scalar perturbations. In this 
context, it is possible to generally parametrize deviations from the 
standard cosmological scenario in the dynamics of perturbations 
by mean of two functions of time and scale introduced in the set 
of Einstein and Boltzmann equations for metric and matter pertur¬ 
bations. More precisely, one can write the Poisson and anisotropy 
equations as follows: 

2 CL^ 

k 'i'==JsUpia,k) , ( 1 ) 

where pA = p6 + 3^(p + P)vis the comoving density perturba- 
tion of matter fields and we have selected the conformal Newtonian 
gauge with and T> representing the perturbation to respectively 
the time-time and space-space diagonal component of the metric. 
And then combine them with the unmodified Boltzmann equations 
for matter fields. 

We shall focus on scalar-tensor theories where the metric 
and the additional scalar degree of freedom obey second order 
eq uations of motion and will a dopt the parametrization introduced 
in teertschinger & Zukinl2008h (BZ) to describe the corresponding 
form of i.e.: 


_ 1 + /3i Al k‘^a^ 

^ 1 + Af 

_ 1 + 1 ^ 2 X 2 

^ 1 +AiFo'* 


( 2 ) 


^ http://www.sfu.ca/~aha25/MGCAMB.html 
^ http://camb.info 
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where we have adopted the convention of Jzhao et alJl2009h and 
^ 1,^2 are dimensionless contants representing couplings, Ai,A 2 
have dimensions of length and s > 0 to ensure that at early times 
GR is recovered. This parametrization gives a very good represen¬ 
tation of scalar-tensor theories in the quasi-static regime, where 
time derivates of the perturbations to the metric and scalar degree of 


freedom are ne glected with respect to their spatial gradients on sub- 
horizon scales dzhao et al .l2009l:lHo iiati et al. l201 j : lde Felice et al.l 


I 2 OIQI : lAmendola et ahll^lSl : ISilvestri et al. 2013h . This is a good 

approximation given the observables that we are considering. Ad¬ 
ditionally, (O sets the evolution of the characteristic lenghtscales 
of the models to a power law in the scale factor. This is of course 
a choice of parametrization for the time dependence of the mass 
scale of the scalar degree of freedom, and other choices are possi¬ 
ble. Nevertheless, as we will discuss in the following, it is a good 
approximation for several scalar-tensor models, and data are not 
that sensitive to the specific choice of the time dependence. 

Equations ^ are built-in in MGCAMB and allow to easily ex¬ 
tract predictions for scalar-tensor models on a ACDM background 
for different observables, including the growth rate. 


2.2 Chameleon models 

Chameleon models are a class of scalar-tensor theories for which 
the additional scalar field has a standard kinetic term and is confor¬ 
mally coupled to matter fields as follows: 

- v(0) 

+Si , (3) 

where is the coupling between the scalar field (j) and the i- 

th matter species. The coupling(s) in general can be a non-linear 
function(s) of the field 0; however, since the value of the field 0 
typically does not change significantly on the time scales associ¬ 
ated to the epoch of structure formation, we will assume it to be 
linear in 0. Since we are dealing with clustering of matter in the 
late universe, it is safe to consider one coupling, i.e. to dark matter; 
that amounts to neglecting differences between baryons and dark 
matter, or simply neglecting baryons, which is safe for the observ¬ 
ables under consideration. 

In the quasi-static regime, (/x, 7 sHp) for Chameleon-type theo¬ 
ries can be well represented by a simplified version of 0 for which: 



Chameleon models under the BZ parametrization, such theoreti¬ 
cal prior has not been gener ally imposed and a wider range of /3i 
has b een explored (see e.g. (iHoiiati et al .1120111 : IPi Valentino et al.l 
l2012h '). Hence, in our analysis we will consider both the case 
with /3i > 1 and the case for which /3i is allowed to be smaller 
than unity, to facilitate comparison. We will refer to the former as 
the Chameleon model, and the latter as the extended Chameleon 
model (eChameleon). We shall emphasize that we consider the 
eChameleon as a purely phenomenological model within with¬ 
out linking it to action 0 , since it would not be viable case of 
the latter. While the eChameleon might correspond to a very spe¬ 
cial subcase of the parametrization 0 , it still represents a possible 
choice for (/x, 'ysiip) and, as we will discuss in Section [6l it will be 
interesting to see what data can say about it. 


2.3 f{R) gravity 

f{R) theories of gravity correspond to the simple modification of 
the Einstein-Hilbert action by the addition of a nonlinear func¬ 
tion of the Ricci scalar. In the past decade they have been exten- 
sivel y explored as candidate m o dels for cosmic accek ration (see 
e.g. jSilvestri & TroddenI l2009l : Ide Felice et al.l l201Clh and refer¬ 
ences therein). They represent a subcase of the larger class of mod¬ 
els described by action 0, corresponding to a universal fixed cou¬ 
pling ai = 7 / 2/3 (j) and are therefore well represented in the quasi¬ 
static regime by the functions JJ]) and conditions (|4]). However, the 
fixed coupling ai = implies that /3i = 4/3 and viable 

f{R) models that closely mimic ACDM have been sh own to cor¬ 
respond to s ~ 4 dZhao et al.ll2009l : lHofiati et al.ll2012h . Therefore 
the number of free parameters in Eqs. JJ]) can be effectively reduced 
to Al, which is then expressed in terms of Bq. The latter is in fact 
the only free parameter needed to label the family of f{R) models 
that reproduce a given expansion history, in our case the ACDM 
one, and can be usuafiy re c onstructed via the so-calle d designer 
approach dSong et al.ll2007l : IPogosian & Silvestrill2008h . Alterna¬ 
tively, one could adopt the recently developed EFTCAMB package 
for an exact implementation of designe r f(R) models that does not 
rely o n the quasi-static approximation dHu et al.ll20l4lRaveri et al.l 
I 2 OI 4 I) [j. The latter method allows to choose different bacgkround 
histories, however for the data and cosmology involved in our anal¬ 
ysis, MGCAMB provides enough accuracy. 


2.4 Growth index parametrization of the growth rate 




2 


\2 


1 ^ s ^ 4 


(4) 


Therefore the effects of Chameleon-type theories on the dy¬ 
namics of linear scalar perturbations on sub-horizon regimes can 
be described with good accuracy in terms of three parameters: 
{/3i, Al, s}. The last condition in © is br oadly valid for m odels 
with runaway and tracking type potentials dZhao et al.ll2009l) . Fol¬ 
lowing a convention which is commonly used for f{R) theories, 
let us express the lengthscale Xf in terms of a new parameter Bq, 
which corresponds to the value of the inve rse mass scale today in 
units of the horizon scale dSong et al.ll2007h : 


Bo = 




(5) 


so that we will work with {f3i, Bo, s}. 

Let us notice that Chameleon theories as defined in action 
13 ), have necessarily /3i ^ 1. However, in previous analysis of 


In the cosmological concordance model, as well as in non¬ 
clustering dark energy models, the growth rate of structure is well 
approximated by: 

a In a 

where Ctmia) = pm{cL)/3MpH^(a), pm is the background 
density of matter and 6 m = 6 pm I Pm- This inspired the fol- 
lowi ng parametrization for deviations in the growth of struc- 

ture dWang & Steinhardtil 19981 : lLindeill2005l : lLinder & Cahnll200% 

f^^miay (7) 

where 7 is commonly referred to as growth index (not to be con¬ 
fused with the 7s Hp defined above, which represents instead the 
gravitational slip). 

^ http://wwwhome.lorentz.leidenuniv.nl/~hu/codes/ 
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Figure 1. The measured fag from different surveys covering redshift range 
0.06 < 2 ; < 0.8. The empty markers represent the reported measurement 
of fag and the hlled markers are for the corrected values for Planck Com- 
sology. The red band shows the Planck ACDM Icr prediction. 

The idea behind this parametrization is that of capturing in¬ 
dependently in Qm and 7 the information from, respectively, the 
expansion and the growth history. Since in our analysis we fix the 
background to ACDM, Qrn(a) is determined by that and the only 
parameter of interest will be 7. While for models of modified grav¬ 
ity and clustering dark energy in general 7 will be a function of 
time and scale, in several cases for the regime of interest it can still 
be safely approximated by a constant, which can di ffer significantly 
from the ACDM value. See (iLinder & Cahnll200% for more details 
and some forms of 7 in alternative theories of gravity. 

In our analysis we will assume 7 is constant and explore con¬ 
straints on it after extracting predictions for the CMB and growth 
of structure from MGCAMB. 


Table 1. Measurement of f(z)ag{z) from various galaxy redshift surveys 
covering redshift between 0.06 to 0.8. 


z 

fcrsiz) 

l/k[h /Mpc] 

Survey 

0.067 

0.42 ± 0.05 

16.0-30 

6dFGRS(2012) 

0.17 

0.51 ±0.06 

6.7-50 

2dFGRS(2004) 

0.22 

0.42 ± 0.07 

3.3-50 

WiggleZ(2011) 

0.25 

0.35 ± 0.06 

30 - 200 

SDSS LRG (2011) 

0.37 

0.46 ± 0.04 

30 - 200 

SDSS LRG(2011) 

0.41 

0.45 ± 0.04 

3.3-50 

WiggleZ(2011) 

0.57 

0.462 ± 0.041 

25 -130 

BOSS CMASS 

0.6 

0.43 ± 0.04 

3.3-50 

WiggleZ(2011) 

0.78 

0.38 ±0.04 

3.3-50 

WiggleZ(2011) 

0.8 

0.47 ±0.08 

6.0-35 

Vipers(2013) 


WMAP7 teennett et al.l2013h likelihood in the analysis. To be able 
to use this fag measurement we need to account for the transfor¬ 
mation to the Planck best fit cosmology (jPlanck Collaboration et al.l 
l2014ah . 

3.2 2dFGRS 

The 2dFGRS (2 degree Field Galaxy Redshift Survey) obtained 
spectra for 221414 galaxies in visible band on the southern sky 
(IColless et al.ll2003h . The survey covers redshift range 0 < z < 
0.25 and has a effective area of 1500 square degree. The RSD mea¬ 
surement was obtained by linearly modeling the observed distortion 
after splitting the ov er-density into radial and angular components 
(jPercival et al.l2004l) . The parameters were fixed at different values 
Us = 1.0 ,iTo = 72. The results were marginalized over power 
spectrum amplitude and bag. We are not using this measurement 
in our analysis for two reasons. First, the survey has a huge over¬ 
lap with 6dFGRS which will lead to a strong correlation between 
the two measurements. Second, the cosmology assumed is quite 
far from WMAP7 and Planck which may cause our linear theory 
approximation used to shift the cosmology might fail. 


3 OBSERVATIONS 

We use meas urements of CMB angular power sp ectrum (Ci) from 
Planck 2013 (jPlanck Collaboration et al.l2014bh combined with the 
measurement of f{z)ag{z) from various redshift surveys covering 
between z = 0.06 to 2 ; = 0.8 listed in Table [T] as our main data 
points. The Figure [T] shows the measurements used with and with¬ 
out corrections and Planck 2013 prediction. We briefiy describe 
each of the surveys and fag measurements in the following sec¬ 
tions. 


3.1 6dFGRS 


The 6dFGRS (6 degree Field Galaxy Redshift Survey has ob¬ 
served 1 25000 galax i es in n ear infrared band across 4/5th of south¬ 
ern sky Ijones et ^ (I2OO9I) . The surveys covers redshift range 
0 < z < 0.18, and has an effective volume equivalent to 
2dFGRS jPercival et al.ll2004h galaxy survey. The RSD measure¬ 
ment was obtaine d using a subsampl e of the survey consisting 
of 81971 galaxies (iBeutler et al.ll201^ . The measurement of fag 
was obtained by fitting 2D correlation function using streaming 
model and fitting range 16- 30 Mpc/h. The Alcock-Paczynski effect 
lAlcock & Paczvnso ( 1979h has been taken into account and it has a 
negligible effect Jb cutler et al.ll2012h . The final measurement uses 


3.3 WiggleZ 

The WiggleZ Dark Energy Survey is a large scale galaxy red¬ 
shift survey of bright emission line galaxies. It has obtained spec¬ 
tra for nearly 200,000 galaxies. The survey covers redshift range 
0.2 < z < 1 . 0 , covering effective area of 800 square degrees 
of equatorial sky (iBlake et al.ll2011bh . The RSD measurement was 
obtained using a sub-sample of the survey consisting of 152,117 
galaxies. T he final result was ob tained by fitting the power spec¬ 
trum using Ijennings et al.l (l 201 lh model in four non-overlapping 
slices of redshift. The measured growth rate is fag{z) = (0.42 zb 
0.07, 0.45 zb 0.04, 0.43 zb 0.04, 0.38 zb 0.04) at effective redshift 
z = (0.22, 0.41, 0.6, 0.78) with non-overlaping redshift slices of 
Zsiice = ([0-1 5 0-3], [0.3, 0.5], [0.5, 0.7], [0.7, 0.9]) respectively. 
We can assume the covariance between the different measurements 
to be zero because they have no volume overlap. 


3.4 SDSS-LRG 

The Sloan Digital Sky Survey (SDSS) data release 7 (DR7) is 
a large- scale galaxy redshift survey of Luminous Red Galaxies 
(LRG) teisenstein et al.ll201lh . The DR7 has obtained spectra of 
106,341 LRGs, covering 10,000 square degree in redshift range 
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0.16 < z < 0.44. The RSD measurement was obtained by mod¬ 
eling monopole and quadruple moment of galaxy auto-correlation 
function using linear theory. The data was divided in two redshift 
bins: 0.16 < z < 0.32 and 0.32 < z < 0.44. The measurements 
of growth rate are f(Js{z) = (0.3512 ± 0.0583, 0. 4602 zb 0.0378) 
at effe ctive redshift of 0.25 and 0.37 respectively jSamushia et al.l 
l2012ah . These measurements are independent because there is no 
overlapping volume between the two redshift slices. 


3.5 BOSSCMASS 


Sloan Digital Sky Su rvey (SPSS) B ary on Oscillation Spectro¬ 
scopic Survey (BOSS: IPa-Wson et ^ (l2013h ) targets high redshift 
(0.4 < z < 0.7) galaxies using a set of color-magnitude cuts. The 
growth rate measureme nt uses the CMASS (Reid et. al. in prep, 
lAnderson et al. ( ^14ah ) sample of galaxies from Data Release 11 
(lAlam et ~ 2015al) . The CMASS sample has 690,826 Luminous 
Red Galaxies (LRGs) covering 8498 square degrees in the redshift 
range 0.43 < z < 0.70, which correspond to an effective vol¬ 
ume of 6 Gpc^. The fag is measured by modeling the monopole 
and quadruple moment of galaxy auto-corre lation using Convolu - 
tion Lagrangian Perturbation Theory fCLPT: [Carlson et aL ( 2013 )) 
in combination with Gaussian Streaming model Jwang et al' 2014 ). 
The reported measurement of growth rate is fa g = 0.462 zb 0.041 
at effective redshift of 0.57 (lAlam et al.l2015bh . 

We are also using the combined measurement of growth rate 
(fcrg), angular diameter distance (Da) and Hubble constant (H) 
measured from the galaxy auto correlation in CMASS sample at 
an effective redshift of 0.57 (lAlam et al .l2015hh . The measurement 
and its covariance are given below and it’s called eCMASS. 


VIPERS 
WiggleZ 
WiggleZ 
CMASS 
> WiggleZ 
SDSS LRG 
SDSS LRG 
WiggleZ 
6dFGRS 


6dFG 

Wigg|SDSS 

5D55 

Wigg 

CMASWigg 

Wigg 


















2.1% 

1.5®/ci 






7.3°J 


1.5%i 













1.2m 



















0.067 0.22 0.25 0.37 0.41 0.57 0.6 0.78 0.8 

2 ; 


Figure 2. Correlation matrix between all the measurements used in our 
analysis. We have estimated the correlation as the fraction of overlap vol¬ 
ume between two survey to the total volume of the two surveys combined. 


CMB anisotropy which is used to measure the cosmic variance 
limited angular power spectrum Cj^ at the last scattering sur¬ 
face. The Planck measurements helps u s constrain the background 
cosmo lo gy to u nprecedented precision JPlanck Collaboration et al.l 
l2014allcll2015al) . We are using the CMB measurements from Planck 
satellite in order to constrain cosmology. We have assumed that 
Planck measurements is independent of the measurement of growth 
rate from various galaxy redshift surveys. 


fag = 0.46, Da = 1401, H = 89.15 


CeCMASS = 


f 0.0018 -0.6752-0.126n 
-0.6752 550.61 45.881 
v-0.1261 45.881 14.019 


( 8 ) 


3.6 VIPERS 

VIMOS Public Extragal actic Redshift Survey 
VIPERS Ide la Torre et al.l (l2013l) ) is a high redshift small 
area galaxy redshift survey. It has obtained spectra for 55,358 
galaxies covering 24 square degree in the sky from redshift range 
0.4 < z < 1.2. The measurement of growth factor uses 45,871 
galaxies covering the redshift range 0.7 < z,1.2. The fag 
measurement is obtained by modeling the monopole and quadruple 
moments of galaxy auto-correlation function between the scale 6 
h“^Mpc and 35 h“^Mpc. They have reported fag = 0.47 zb 0.08 
at effective redshift of 0.8. The perturbation theory used in the anal¬ 
ysis has been tested against W-body simulation a nd shown to work 
at mi ldly non-linear scale below 10 h“^Mpc dde la Torre et al.l 
l2013h . 


3.7 Planck CMB 

Planck is a space mission dedicated to the measurement of CMB 
anisotropies. It is the third-generation of all sky CMB experiment 
following COBE and WMAP. The primary aim of the mission is 
to measure the temperature and polarization anisotropies over the 
entire sky. The Planck mission provides a high resolution map of 


3.8 Correlation Matrix 

We use the measurements of fag from 6 different surveys. Al¬ 
though these surveys are largely independent, and in some cases 
they probe different biased tracers, they are measuring inherently 
the same matter density field. Therefore, the parts of the survey ob¬ 
serving same volume of sky cannot be treated as independent. We 
have predicted an upper limit to the overlap volume using the data 
from different surveys. We have estimated the fractional overlap 
volume between any two samples as the ratio of the overlap vol¬ 
ume to the total volume of the two samples. We estimate the corre¬ 
lation between two measurements as the fractional overlap volume 
between the two measurements. Pigure|2] shows our estimate of the 
correlation between the surveys. The four measurements of Wig¬ 
gleZ survey cover the redshift range between 0.1 to 0.9 and hence 
show most correlation with other measurements like SDSS LRG 
and CMASS in the same redshift range. 


4 POTENTIAL SYSTEMATICS 

The collection of fag data points that we are using in this analy¬ 
sis contain measurements from several different surveys, obtained 
during the last decade, each with a different pipeline. Eurthermore, 
often the latter implicitly assumes a GR modeling, which does not 
take into account the different predictions for the growth factor in 
modified theories of gravity. It is important to account for some cru¬ 
cial differences in order to use these measurements in our analysis. 
We have looked at following different aspects of measurements and 
theoretical prediction before using them in our analysis. 
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4.1 Fiducial Cosmology of the growth rate (/erg) 


The measurements of fag have been obtained over the 
time when we had trans ition from WMAP best fit cosmol- 
Qgy (iHinshaw et aP 20131) to th e Planck best fit cosmology 
(jPlanck Collaboration et alJ 2014ah. Since w e are using Planck 
likelihood dPlanck Collaboration et al.ll2014bh in our analysis, we 
have decided to convert all the measurements to Planck cosmology. 
The 3 dimensional correlation function can be transformed from 
WMAP to the Planck cosm ology using Alcock-Paezynski effect 
(lAlcock & PaczvnsO (Il979h ). 


^pianck{r \\, r±, (j)) = ^wMAp{a\\r \\, 0), (9) 

where is the ratio of the bubble parameters (ay = 
HwMAPIHpianck) and a± is the ratio of the angular diameter 
distances ( = jyp^rick jjjw map sepa¬ 

rations along the line of sight and perpendicular to the line of sight 
and 0 is the angular position of pair separation vector in the plane 
perpendicular to the line of sight from a reference direction. In prac¬ 
tice the correlation function is isotropic along 0. We can calculate 
the corresponding power spectrum by applying Fourier transform 
to correlation function. 


f’p(ancfe(fc||,fcx,fc<A) = J driidr±r±d(l>^pianck{r\\,r±,(l>)e 

( 10 ) 

= f dr\\dr'j_ — ^^CwMAp{r\\,r'j_, (j))e ' 

( 11 ) 

_ PwMAp{k\\/a\\,k±/a±,k^) 

2 

The Kaiser formula for RSD gives the redshift space correla - 
tion function as Pg(k,fi) = b‘^Pm(k)(l + (lKaisejll98% . 

Using the linear theory kaiser prediction and the above approxi¬ 
mation between WMAP and planck power spectrum, we can get a 
relation to transform the growth function from WMAP to planck 
cosmology. 


1 fiplanckl^j ^ 

1 + PwMAPp ^ 

/ Pplanck{P ■) /^^) 

' PwMAp{k, /jf 

(13) 

= C7. 

l^-T 

(14) 

\ 

1 anal 



where C is the ratio of isotropic matter power spectrum with 
WMAP and planck cosmology integrated over scale used in [3 mea¬ 
surement. 





WMAP 


■ planck 


(fc') 


(15) 


When right hand side of eauaion dldl) is close to 1, then we can 
approximate the above equation as follows: 


^planck = ^WMAPC-^^ ^ 


(16) 


2 

The ratio can be obtained using simple trigonometry 

/i, 

which gives following equations, where the last equation is aprrox- 
imation for a y ^ a\. 



We can substitute equation (fTTl) in equation ([T^ in order to 
get the required scaling for /(growth factor) assuming that bias 
measured is proportional to the ag of the cosmology used. 


^planck 


PWMAPC 



(18) 


/ ay \ / ^planck \ ^ 

f^^planck f^^WMAP^ y J ( ^-W M AP 1 (19) 

We have tested prediction of equ ation ilSL against th e mea¬ 
surement of fag reported in Table 2 of I Alam et al.l (l2015bh at red- 
shift 0.57 using both Planck and WMAP cosmology. 


4.2 Scale dependence 

General Relativity predicts a scale independent growth factor. One 
of the important features of the modified gravity theories we are 
considering is that they predict a scale dependent growth factor 
which has a transition from high to low growth at certain scale 
which depends on the redshift z and the model parameters. The 
mesurements we use from the different surveys, assumes a scale- 
independent fag and uses characteristic length scale while analyz¬ 
ing data. In order to account for all these effects we have done our 
analysis in two different ways. In the first method, we assume that 
the measurements corresponds to an effective k and in the second 
method, we treat the average theoretical prediction over range of k 
used in fag analysis. 

Figures [TT] and [12] show the parameters constraint for 
Chameleon models and f{R) gravity. The grey and red contours 
result from using two different model predictions to test the scale 
dependence. The grey contours correspond to the model where we 
average fag over k used in respective fag analysis and red con¬ 
tours correspond to fag evaluated at /c = 0.2h Mpc“^. It is evi¬ 
dent from the plots that, at the current level of uncertainty, we ob¬ 
tain very similar constraint and hence do not detect any significant 
effect of scale dependence of fag . 


4.3 Other systematics 

The measurements of fag are reported at the mean redshift of the 
surveys. But the galaxies used have a redshift distribution which 
in principle can be taken into account by integrating the theoret¬ 
ical prediction. This should be a very small effec t because the 
fag is relatively smooth and fiat (see Figure ?? and iHuterer et al.l 
(l2013h )for the redshift range of the survey and also because the 
survey window for every individual measurement is small. An¬ 
other important point is the assumption of GR based modeling for 
the measurement. We have looked at the modeling assumption for 
each of the measurements. All measurements of fag except Wig- 
gleZ and VIPERS , allow the deviation from GR through Alcock- 
Paezynski effect (lAlcock & PaczvnsO Il979h which justifies our 
use of modified gravity models. The inclusion of AP in WiggleZ 
and VIPERS will marginally increase the error on the measure¬ 
ments. Different surveys use different ranges of scale in the RSD 
analysis. This will be important especially while analyzing mod¬ 
ified gravity models. To account for the different scales used we 
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evaluate the prediction for each survey averaged over the scale used Table 2. The list of extension parameters for all the models used in our 

in the respective analysis. analysis. For each parameter we provide their symbol, prior range, central 

value and Icr error. 


5 ANALYSIS 

We have measurement of fag from various surveys covering red- 
shift range 0.06-0.8 (see Table: [T])- We first correct these measure¬ 
ments for the shift from WMAP cosmology to planck cosmology 
as described in section dd.ll) . The next step is to evaluate prediction 
from different modified gravity theories by evolving a full set of 
linear perturbation equations. The theoretical predictions for fag 
is generally scale and redshift dependent (see section There¬ 
fore, we consider two cases for theoretical prediction: 1) evaluate 
fag at effective k and 2) evaluate fag averaged over range of k 
used in measurements. We also predict for different modified 
gravity theories. Finally we define our likelihood, which consists of 
three parts one by matching planck temperature fiuctuation , 
second by matching growth factor from Table [T] and third by us¬ 
ing eCMASS data as shown in equationfS] Therefore, we define the 
likelihood as follows: 

JC = CplanckjCfagCeCMASS ( 20 ) 

( 21 ) 

X/<T8 = A/o-sC'^A/ctI’ (22) 

The A fag is the deviation of the theoretical prediction 
from the mesurement and C~^ is the inverse of covariance 
which has diagonal error for different surveys and correlation be¬ 
tween measurement as described in section B.SI . Note that we 
do not include fag from CMASS while using eCMASS with 
fag{z) to avoid double counting. This likelihood is sampled 
using modifi ed version of COSMOMC (iLewis & Bridlel 120021 : 
iHoiiati etaP l201lh . We sample over 6 cosmology parameters 
, 1OO0MC, T, risJog (10^ ^ As)} and all 18 planck 
nuisanc e parameters as described in IPlanck Collaboration et al.l 
(l2014ah with the respective extension parameters or modified grav¬ 
ity parameters. The priors we have used on all the parameter s are 
the same as the priors in IPlanck Collaboration et al.l (l2014ah and 
the priors we used on the parameters of modified gravity model are 
given in Table m 


Model 

Parameter 

prior range 

posterior 

mCDM 

mo 

-2.0 - 0.0 

-0.873 ± 0.077 

WQWaCDM 

mo 

-2.0 - 0.0 

-0.943 ±0.168 


Wa 

q 

q 

0.156 ±0.361 

oACDM 


b 

b 

1 

-0.0024 ± 0.0032 


/3i 

0-2.0 

1.23 ±0.29 


/32 

0-2.0 

0.93 ± 0.44 

Scalar-tensor 

X 10“® 

0-1 

0.49 ± 0.29 


A| X 10-® 

0-1 

0.41 ±0.28 


s 

q 

q 

2.80 ±0.84 


I3i 

1.0-2.0 

< 1.008 

Chameleon 

Bo 

q 

o 

< 1.0 


s 

q 

q 

2.27 < s < 4.0 


I3i 

0-2.0 

0.932 ±0.031 

eChameleon 

Bo 

0- 1.0 

< 0.613 


s 

q 

q 

2.69 < s < 4.0 

f(R) 

Bo 

10“'^®’- 10-"^ 

< 1.32 X 10^ 

Growth index 

7 

0.2 - 0.8 

0.611 ±0.072 



6 RESULTS 

We have combined CMB data set and measurements of 
growth from various redshift surveys in order to constrain 
the parameters of standard cosmology (ACDM), extended 
cosmology models and modified gravity. Our analysis gives 
consistent constraints for the standard ACDM parameters 
{Qbh^ , 1OO0MC, T, ris, log {10^^ As)} as shown inTable[3] 
Figure|4] shows the constraint on flrn-cfg plane for ACDM, wCDM, 
oACDM, Scalar-tensor model. Chameleon gravity, eChameleon, 
f{R) and growth index parametrization. Theses are our best con¬ 
straints obtained using planck -i-eCMASS -i- fag{z). 


6.1 ACDM 

Figure [5] shows the one dimensional marginalized likelihood for 
standard ACDM cosmology. The black line shows the contraints 


Figure 3. The black points show the corrected fag used in our analysis, 
along with the errorbar. Lines of different colors show the best fit for the 
various models used in our analysis. The best fit and are for the case of 
Planck +fag + eCMASS fits. Notice that the eChameleon model predicts 
the smallest growth rate by preferring lower values of the coupling constant 
(/3i), even though the scalar amplitude of primordial power spectrum is 
high. 


from Planck 2013 alone. The red ,blue and magenta lines are pos¬ 
terior obtained for the data set combinations planck-i-eCMASS, 
p\sLnck+fag{z) and p\sLnck+QCMASS+fag{z) respectively. Our 
parameter constraints are completely consistent with the Planck 
2013 results. Adding measurements of the growth rate to Planck 
data does not improve the results (see Figure |5]) due to already tight 
constraints from Planck observations (see Figure [T]). 
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[bp!] 


Table 3. The list of standard ACDM parameters used in our analysis. For each parameter we provide its symbol, prior range, central value and la error. We 
have used the same prior as Planck2013 on these parameters. We have also marginalized over all the nuisance parameters of Planck likelihood. We report the 
results for each of the model analyzed in this paper. 

Models 1006»mc t ns In(lO^OAs) 


prior range 

ACDM 

mCDM 

moiCaCDM 

oACDM 

Scalar-tensor 

Chameleon 

eChameleon 

f{R) 

growth index (7) 


0.005-0.10 
0.0219 ± 0.0002 
0.0221 ± 0.0003 
0.0221 ± 0.0003 
0.0220 ± 0.0003 
0.0221 ± 0.0003 
0.0219 ± 0.0002 
0.0218 ± 0.0003 
0.0221 ± 0.0004 
0.0218 ± 0.0003 


0.001-0.99 
0.1208 ±0.0020 
0.1183 ±0.0028 
0.1181 ±0.0029 
0.1191 ±0.0031 
0.1199 ±0.0020 
0.1205 ±0.0020 
0.1222 ±0.0023 
0.1182 ±0.0033 
0.1214 ±0.0023 


0.50-10.0 
1.0410 ±0.0006 
1.0414 ±0.0006 
1.0415 ±0.0007 
1.0413 ±0.0007 
1.0412 ±0.0006 
1.0411 ±0.0006 
1.0409 ±0.0006 
1.0414 ±0.0008 
1.0409 ±0.0006 


0 . 01 - 0.8 

0.0442 ± 0.0236 
0.0911 ±0.0449 
0.0906 ± 0.0454 
0.0518 ± 0.0278 
0.0333 ± 0.0198 
0.0390 ± 0.0222 
0.1313 ±0.0467 
0.0733 ± 0.0354 
0.0699 ± 0.0400 


0.9-1.1 

0.953 ± 0.0068 
0.9615 ± 0.0097 
0.9619 ± 0.0099 
0.9582 ± 0.0094 
0.9591 ± 0.0071 
0.9539 ± 0.0067 
0.9537 ± 0.0079 
0.9607 ±0.0101 
0.9525 ± 0.0075 


2.7-4.0 

3.0007 ± 0.0450 
3.0884 ± 0.0843 
3.0871 ± 0.0850 
3.0118 ±0.0514 
2.9769 ± 0.0377 
2.9894 ± 0.0425 
3.1780 ±0.0914 
3.0526 ± 0.0663 
3.0534 ± 0.0788 



Figure 4. The hgure shows the Icr and 2f7 regions for each of the models 
considered in this paper in plane. It shows that the posterior like¬ 

lihood is consistent for each of the model in this parameter space. The top 
plots shows the models which are extension to ACDM and the bottom plot 
shows the modified gravity models. 





Figure 5. ACDM: We use GR as the model for gravity to determine the 
growth factor and ht for fagiz) and eCMASS measurement with Planck 
likelihood. The black line shows the contraints from Planck 2013 alone. The 
red ,blue and magenta lines are posterior obtained for the data set combi¬ 
nations planck-i-eCMASS, plmck+fag (z) and planck+eCMASS+/(78 (z) 
respectively. The two most prominent effect are in optical depth r and scalar 
amplitude of primordial power spectrum As . Which is also reflected in the 
derived parameter ag and mid redshift of re-ionization Zre- 

6.2 Dark Energy Equation of state (ti;CDM) 

We have looked at the it;CDM, i.e the one parameter extension of 
ACDM where the dark energy equation of state is a eonstant, w. 
Figure [6] shows the two dimensional likelihood of wo and Qm- The 
grey eontours are Planck only constraint (wo = —1.27 ± 0.42), 
red contours are Planck and eCMASS (wo = —0.92 ± 0.10) and 
blue contours show Planck combined with eCMASS and growth 
factor measurements (wo = —0.87 ± 0.077). We obtain wo = 
—0.87 ± 0.077 (8.8% measurement) which is consistent with the 
fiducial value of u; = — 1 for ACDM. The constraint we obtained 
is similar in precision as compared to BAO only, but has different 
degeneracy. Therefore combined measurement of growth rate and 
anisotropic BAO for all of these surveys will help us improve the 
precision of wo. 

6.3 Time-dependent Dark Energy (i(;otUaCDM) 

The tuCDM model which proposes a constant dark energy is lim¬ 
ited in its physical characteristics. Many models propose time- 
dependent dark energy which is popularly tested using linear re¬ 
lation w(z) = Wo + with Wo and Wa as free parameters. 

This model has been shown to match exact so lutions of distance, 
Hubble, growth to the 10“^ level of accuracy jde Putter & LindeJ 
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Figure 6. wCDM: The two dimensional posterior likelihood w and Q-m 
for mCDM. The grey contour is for Planck (wq = —1.27 ± 0.42), red 
contour is combined constraint from Planck and eCMASS (mo = —0.92 zb 
0.10). The blue contour represents constraint from combining Planck with 
eCMASS and fagiz) (wq = -0.87 ± 0.077). 



Figure 8. oACDM: The two dimensional posterior likelihood of Qk and 
Qrn for oACDM. The grey contour is for Planck = —0.060 ± 0.047), 
red contour is combined constraint from Planck and eCMASS = 
—0.0024 zb 0.0034). The blue contour represents results from combining 
Planck with eCMASS and fcrgiz) (flk = -0.0024 ± 0.0032). 



Figure 7. wOwaCDM: The two dimensional posterior likehhood of mo 
and Wa for time-dependent dark energy model. The grey contour is for 
Planck (mo = —0.99 ± 0.52, Wa = —1.50 ± 1.46), red contour is com¬ 
bined constraint from Planck and eCMASS (mo = —1.23 ± 0.26, Wa = 
0.63 ± 0.49). The blue contour represents results from combining Planck 
with eCMASS and fcrgiz) (mo = -0.94 ± 0.17, Wa = 0.16 ± 0.36). 

I 2 OO 8 I) for a wide variety of scalar field (and modified gravity) mod¬ 
els. The dynamical evolution of w{z) can change the growth factor 
significantly and leave an imprint on the CMB. The combination 
of CMB and collection of growth factor at different redshifts is a 
unique way to test the time-dependent dark energy model. 

Figure |7] shows the Icr and 2cr region for (mo, Wa). The grey 
contour is from the Planck temperature power spectrum data alone 
(mo = —0.99 zb 0.52, Wa = —1.50 zb 1.46). The red contours are 
from Planck and eCMASS (mo = —1.23 zb 0.26, Wa — 0.63 zb 
0.49) and blue contour shows planck combined with eCMASS and 


growth factor measurement (mo = —0.94 zb 0.17, ma = 0.16 zb 
0.36). The ACDM prediction of (mo, ma) = (—1,0) is completely 
consistent with our posterior. We have obtained constraint on mo = 
—0.94 zb 0.17 (18% measurement) and 1 + ma = 1.16 zb 0.36 
(31% measurement) which is stronger con straint than the curren t 
best measurement of ma = —0.2 zb 0.4 from lXubourg et ^(l2014l) . 

6.4 Spatial Curvature (oACDM) 

We consider a model with spatial curvature parametrized with VLk 
as free parameter called oACDM along with ACDM parameters. 
Figure [8] shows the posterior for the Q.k and Qm plane. The grey 
contour is from the Planck temperature power spectrum data alone 
(Qfc = —0.060 zb 0.047). The red contours are from Planck and 
eCMASS (Qfc = —0.0024 zb 0.0034) and blue contour shows 
Planck combined with eCMASS and growth factor measurements 
(Q^ = —0.0024 zb 0.0032). The ACDM prediction of Qfc = 0 
is completely consistent with our posterior. We have obtained con¬ 
straint on 1 + = 0.9976 zb 0.0032 (0.3% meas urement) which 

is com petitive with the current best measurements (lAubourg et al.l 
I 2 OI 4 I) . It will be interesting to see if combined RSD and BAO at all 
redshifts will give any improvement on the precision of curvature. 

6.5 Scalar-Tensor Gravity (BZ parametrization) 

The general Scalar-tensor theories of gravity is analyzed using five 
parameter model called BZ parametrization . The five parameters 
of Scalar-tensor gravity (/3i,/32, Ai, A 2 , s) are constrained along 
with the standard ACDM parameters using Planck, fag (z) and eC¬ 
MASS measurements. BZ model predicts a scale dependent growth 
rate (fag{k, z)), whereas the measurements are at some effective 
k. In order to incorporate the k-dependence in our analysis, we 
use the two different approaches described in Section 14.21 Fig¬ 
ure |9] shows two dimensional posterior in the plane (I3i, ^ 2 ) • The 
green contour is combined constraint from Planck and eCMASS 
(f3i = 1.18 zb 0.29,/32 = 0.95 zb 0.43). The grey contour is 
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Figure 9. BZ; The two dimensional posterior likelihood of ^ 1-^2 for five 
parameter Scalar-tensor theory parametrized through the BZ form of EJ. 
The green contour is the combined constraint from Planck and eCMASS 
(13 1 = 1.18 ± 0.29, /32 = 0.95 ± 0.43).The grey contour is the combined 
constraint from Planck and fcrgiz) with averaged over k (/3i = 1.24 ± 
0.3, ^2 = 0.96±0.45), red contour is the combined constraint from Planck 
mdfasiz) at effective k = 0.2/iMpc“^ (^1 = 1.24±0.3,/32 = 0.95± 
0.45). The blue contour represents results from the combination of Planck, 
fas iz) and eCMASS (/3i = 1.23 ± 0.29, ^2 = 0.93 ± 0.44). 



2.76 


2.30r 


1.84 


1.38 


0.92r 


Planck-1- eCMASS at k=0.2hMpc ^ 

Planck -I- fag{z) mean k 
Planck-P fag{z) at k=0.2hMpc^^ 

Planck-l- eCMASS +fag(z) at k=0.2hMpc~^ 


0.46 


0.00 


0.53 


the combined constraint from Planck and fagiz) with averaged 
over k (/3i = 1.24 zb 0.3, p 2 = 0.96 zb 0.45); the red contour 
is the combined constraint from Planck and fasiz) at effective 
k = 0 . 2/1 Mpc"^ (/3i = 1.24 ± 0.3 ,/32 = 0.95 ± 0.45). The 
blue contour represents results from the combination of Planck, 
fasiz) and eCMASS (/3i = 1.23 zb 0.29 ,/32 = 0.93 zb 0.44). 
We obtain the following joint constraint on the five BZ parameters: 
/3i = 1.23 ± 0.29 ,/32 = 0.93 ±0.44, Af(xl0“®) = 0.49 ± 0.29 
,Ai(xl0“®) = 0.41 ± 0.28 and s = 2.80 ± 0.84. By looking 
at the joint 2D likelihood for (/ 3 i,/ 32 ) in Fig. we notice that 
there is a strong degeneracy between the two parameters which re- 
fiects the degeneracy between /x and 7 for the observables that we 
are using. Similar results have be en found in (iHoiiati et al.ll2012l : 
IPlanck Collaboration et al.l l2015bh . For the next models that we 
will discuss, /3i and (32 are not independent and this will allow data 
to place more stringent constraints. 

While the constraints on the length scale of the scalar field 
(Ai, A 2 ) and (s) are very broad, the one on the coupling, and 
p 2 , is the first ever constraint obtained on these parameters for gen¬ 
eral Scalar-tensor gravity. The discrepancy in the strength of the 
constraints on the coupling and on the length scale, can be linked 
to the fact that data strongly prefer values of the coupling constants 
close to 1. For such values, the scale and time dependences in (/x, 7) 
becomes less important and therefore are loosely constrained. We 
will encounter this again in the Chameleon and f(R) gravity cases. 


6.6 Chameleon Gravity 

The three parameters of Chameleon gravity (/3i,Bo,s) are con¬ 
strained along with the standard ACDM parameters using Planck, 
fcTsiz) and eCMASS measurements. Chameleon models predict 
a scale dependent growth rate (fa8{k,z)), whereas the measure¬ 
ments are at some effective k. In order to incorporate the k- 



Figure 10. Chameleon Theory: The two dimensional posterior likelihood 
for Chameleon gravity. The green contour is the combined constraint from 
Planck and eCMASS (/3i < 1.013).The grey contour is the combined con¬ 
straint from Planck and fcrgiz) with averaged over k (/3i < 1.010), the 
red contour is the combined constraint from Planck and fagiz) at effective 
k = 0.2h Mpc“^ (^1 < 1.010). The blue contour represents results from 
the combination of Planck, eCMASS and fcrg(z) (/3i < 1.008). 


Planck-l- eCMASS at k=0.2hMpc~^ 

Planck-I-/cr 8 ( 2 ;) mean fc 
Planck-h fag(z) at k=0.2hMpc ^ 

Planck-l- eCMASS +f(7g(z) at k=0.2hMpc~^ 



Figure 11. eChameleon Theory: The two dimensional posterior likelihood 
of /3i and Qm for extended Chameleon gravity. The green contour is the 
combined constraint from Planck and eCMASS (/3i = 0.932 ± 0.04). 
The grey contour is the combined constraint from Planck and fag (z) with 
averaged over k (/3i = 0.940 ± 0.032), red contour is combined constraint 
from Planck and fag(z) at effective k = 0.2h Mpc“^ (/3i = 0.936 ± 
0.032). The blue contour represents results from the combination of Planck, 
fasiz) and eCMASS ((3i = 0.932 ± 0.031). 


dependence in our analysis, we use the two different approaches 
described in Section lOl FigurefTOlshows the two dimensional pos¬ 
terior in the plane (Qrn^Pi), and (s,/3i). The grey and 

red contours show the posteriors from combined data set of Planck 
and growth rate measurements. The red contours are likelihood 
while evaluating the growth rate at an effective k (/3i < 1.010), 
whereas grey contours are for the case when we use an effective 
growth rate, averaged over the scales used in the actual fag mea¬ 
surement (/3i < 1.010). The green contour is combined constraint 
from Planck and eCMASS (/3i < 1.013). Finally, the blue con¬ 
tours show the posterior from combined data of Planck, eCMASS 
and growth rate (/3i < 1.008). We obtain the following joint con¬ 
straint on the three Chameleon parameters: /3i < 1.008, Bq < 1.0 
and 2.27 < s < 4. While the constraints on the length scale of the 
scalar field, Bq and s is very broad, the one on the coupling, /3i, 
is very strong and predicts = 1 to 0.8%, bringing fi to its GR 
value. As we already discussed for the Scalar-tensor case, the dis¬ 
crepancy in the strenght of these constraints is due to the fact that 
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data prefer values of the coupling constant close to 1 , for which the 
time and scale dependence of {jJL^^sUp) become negligible. This 
is even more the case for Chameleon models, where the theoreti¬ 
cal prior forces /3i > 1 , which corresponds to enhanced growth, 
and data consequently require very small values for this coupling, 
pushing fi very close to its GR value. 

We have also looked at the extended Chameleon model where 
we allow /3i to be less than 1 following previous analysis of this 
model. Figure [TT] shows the two dimensional posterior in the plane 
(Qrn,/^i). The red contours are likelihood while evaluating the 
growth rate at an effective k {/3i = 0.940 zb 0.032), whereas 
grey contours are for the case when we use an effective growth 
rate, averaged over the scales used in the actual fag measurement 
(/3i = 0.936 zb 0.032). The green contour is combined constraint 
from Planck and eCMASS = 0.932 zb 0.04). Finally, the blue 
contours show the posterior from combined data of Planck, eC¬ 
MASS and growth rate (/3i = 0.932 zb 0.031). We obtain the 
following joint constraint on the three eChameleon parameters: 
/3i = 0.932 zb 0.031, Bo < 0.613 and 2.69 < s < 4. Like 
in the more general Scalar-tensor case, while the constraints on 
the length scale of the scalar field. Bo and s are very broad, the 
one on the coupling, ^i, is an huge improvement on the previ¬ 
ous constraint of /3i = 1.3 zb 0.25 (1 9.2 % measuremen t) using 
WMAP CMB, SNe and ISW dataset jHoiiati et alJl201lh . Let us 
notice that when we constrain jointly the three eChameleon pa¬ 
rameters, data select a region in the parameter space which cor¬ 
responds to /3 1 < 1, i.e. to suppressed growth. This region ex¬ 
cludes standard Chameleon models, including f{R) theories, for 
which /3i > 1 and the growth is enhanced. After all, as we have 
seen above and will see in the next Section, the same data place 
very stringent constraints on Chameleon and f{R) models, forcing 
them to be very close to ACDM (see Figure [TT] and O. Hence the 
combination of data sets that we employ favor models with a sup¬ 
pressed growth rate, which adopting the BZ parametrization can be 
obtained with /3i < 1 ; a suppressed growth was favor ed also by 
the data set used in (jPlanck Collaboration et al.ll2015hh . although 
in that case the authors employed a time-dependent parametriza¬ 
tion. Theoretically viabl e scalar-tensor mod els with a suppressed 
growth are discussed in (jPerenon et al.ll2015b . where they are ana¬ 
lyzed via a scale-independent parametrization in the effective field 
theory language. 


6.7 f(R) theory 

We consider one parameter (Bo) model of f(R) gravity. The pa¬ 
rameter Bo parameterizes the deviation from ACDM. The model 
approaches GR when Bo is zero. Similar to Chameleon theory, 
f{R) gravity predicts a scale dependent growth rate (fagik^z)). 
Figure [12] shows the two dimensional posterior in Bo and Qm 
plane. The green contour is combined constraint from Planck and 
eCMASS (Bo < 3.43 x 10“^). The grey and red contours show 
posterior from combined data set of Planck and growth rate mea¬ 
surements. The red contours are likelihood while evaluating the 
growth rate at an effective k (Bo < 1.89 x 10 “^) where as 
grey contours are for the case when we use effective growth rate, 
which is averaged over scales used in the actual fag measure¬ 
ments (Bo < 2.77 X 10“^). The blue contours show the pos¬ 
terior from combined data of Planck, eCMASS and growth rate 
(Bo < 1.36 X 10“^). We obtained Bo < 1.36 x 10“^ (la C.L.), 
which is an improvement by a factor of 4 on the most recent con- 
strai nt fro m large scale structure of Bo = 5.7 x 10“^ (Icr C.L.) 
(IXull2015h . Our constraint is competitive with the constraint from 
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Figure 12. f(R) Gravity; The two dimensional posterior likelihood of 
Bo and Qm for f(R) gravity. The green contour is combined constraint 
from Planck and eCMASS (Bo < 3.43 x 10“^). The grey contour is 
the combined constraint from Planck and fag(z) with averaged over k 
(Bo < 2.77 X 10“^), the red contour is the combined constraint from 
Planck and fag(z) at effective k = 0.2h Mpc“^ (Bo < 1.89 x 10“^). 
The blue contour represents results from the combination of Planck, fag (z) 
and eCMASS (Bo < 1.36 x IQ-^) 


solar system tests and clusters (IHu & Sawickill200^ : ISchmidt et al.l 
l2009l> . 


6.8 Growth index ( 7 ) parametrization 

The standard cosmological model, based on GR, predicts a precise 
value for the growth factor in the linear regime, i.e. / = . In 

order to test deviations from GR, we have parameterized the growth 
factor using growth index 7 ( cite Linder & Cahn 2008) as / = 
. The marginalized two dimensional likelihood for flm and 7 is 
shown in Figure [T^ The grey contour is combined constraint from 
Planck and eCMASS (7 = 0.477 zb 0.096). The red contours show 
the constraint obtained using Planck and fag (z) measurement (7 = 
0.595 zb 0.079) and the blue contours are for combined data set of 
Planck with fag(z) and eCMASS (7 = 0.612 zb 0.072). We have 
obtained 7 = 0.612 zb 0.072 (11.7% measurement) completely 
consistent with the general relativity prediction. 


7 DISCUSSION 

We have constrained the parameters of the standard cosmological 
model, ACDM, as well as those of various extensions using the 
current measurements of growth rate between redshift 0.06 and 
0.83 (Figure [4]), eCMASS and Planck 2013. We have been care¬ 
ful with several important details while combining results from 
various surveys and different cosmologies of measurements. We 
have first showed that the standard ACDM parameter space has a 
consistent posterior, independent of the model considered except 
for Chameleon gravity. Next, we focused on each model and an¬ 
alyzed the constraint on its extension parameters. As for the stan¬ 
dard model, ACDM, using the growth factor we do not improve 
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Figure 13. Growth index ( 7 ): The two dimensional posterior likelihood 
of 7 and for growth index parametrization. The grey contour is the 
combined constraint from Planck and eCMASS (7 = 0.477 ± 0.096). 
The red contour is the combined constraint from Planck and fagiz) (7 = 
0.595 ±0.079). The blue contour represents results from the combination 
of Planck, fag iz) and eCMASS (7 = 0.612 ± 0.072) 


constraints on any of its parameters because the growth rate is al¬ 
ready highly constrained with Planck measurement for the standard 
model of cosmology. It is impressive to notice that ACDM, with¬ 
out any extra parameter, is completely consistent with the measure¬ 
ments of fag from very different galaxy types and redshifts. In the 
case of the extension where the dark energy equation of state is 
constant but free to vary, icCDM, we obtain wq = —0.87 ± 0.077 
(8.8% measurement). This is a 3.7 times improvement on the pre¬ 
cision compared to Planck only measurement w = —1.27 ± 0.42 
(33% measurement) an d comparable to the 8% measurement of 
ISamushia et al.l (12012hl) . Our measurement prefers w < — 1 at Icr 
level. We have also noticed that the growth rate and BAO have 
slightly different degeneracy for wCDM. This shows the potential 
to improve the constraint on w by combining the growth rate and 
BAO measurements from a range of galaxy redshift surveys. How¬ 
ever, one difficulty in doing so, is to model the correlation between 
the measurement of growth rate and BAO. 

We also report one of the best measurements on the parameters 
of the model with a time-dependent equation of state, woWaCDM. 
We have measured wq = —0.94 ±0.17 (18% measurement) and 
l±u;a = 1.16±0.36(31% measurement). This represents a sig- 
nificant improvement on Wg compared to all other me asurements 
(jPlanck Collaboration et al.l2014al : lAubourg et al.l2014h . The mea¬ 
surements of fag, H and Da in eCMASS help to constrain ACDM 
parameters, while the evolution of the growth rate over a large red- 
shift range, obtained through measurements of fag{z) at multi¬ 
ple redshifts, improves the constrain on evolving dark energy. This 
hints at the potential of using combined growth rate and anisotropic 
BAO as function of redshifts , when future surveys like eBOSS and 
Euclid dLaureiis et al.ll20llh will provide much stronger growth 
rate and BAO constraints at much higher redshifts. We have also 
looked at the possibility of a non-zero curvature for the universe, 
oACDM, finding l±Ofc = 0.9976 ±0.0032 (0.3 % measurement), 
which is same as the best constraint reported in dSamushia et al.l 
l2012bl) . We notice that the optical depth (r) and amplitude of scalar 


power spectrum (As) is relatively low for oACDM, which pre¬ 
dicts smaller redshift of reionization (zre = 7.20 ± 2.81) but it 
is abo ve the lower hmit ob served through Lymana Forest observa¬ 
tions dBecker et al.l d200lh ). 

We have also looked at some of the popular modifications of 
gravity and found no significant deviations from GR using growth 
rate and Planck 2013 measurement. We have investigated gen¬ 
eral scalar±enso£_theori£^suMer the parametrization introduced 
in dBertschinger & ZukinI l2008l) . constraining the corresponding 
five parameters. We have found constraints on the two coupling 
parameters (/3i = 1.23 ± 0.29, (32 — 0.93 ± 0.44) while the pos¬ 
terior of other three parameters were largely non-constraining. We 
then restricted to the subset of Chameleon theories, for which only 
three parameters are needed. While imposing a theoretical bound 
of (3 1 > 1, we constrained the coupling of Chameleon theories to 
Pi < 1.008 (1 a C.L.), while jointly varying the remaining two 
free parameters that describe the lengthscale of the scalar degree of 
freedom, {Bo, s}. While the latter are loosely constrained by data, 
the constraint on the coupling is quite stringent. We explored also 
an extension of Chameleon models, that we dubbed eChameleon, 
where we let the coupling Pi vary within the range [0, 2]. Also 
in this case, data place a stringent bound on the coupling, while 
loosely constraining {Bo, s}. Interestingly, for this case data select 
aregion where/3i < 1, with the bound/3i = 0.932±0.031; the lat¬ 
ter corresponds to a region of the parameter space for which growth 
is suppressed. This improves significantly over previous analysis, 
e.g . the bound Pi = 1-3 ± 0.25 (19.2 % measurement) obtained 
in dHoiiati et al'll201ll) using WMAP CMB, SNe and ISW dataset. 
This excludes standard Chameleon models, including /(B) the¬ 
ories, for which Pi > 1 and the growth is enhanced. After all, 
the same data place very stringent constraints on the latter mod¬ 
els, forcing them to be very close to ACDM (see Figure [TT] and 
[13. We also notice that the optical depth (r) and amplitude of 
scalar power spectrum (As) is higher for eChameleon gravity. This 
predicts higher redshift of reionization (zre = 14.43 ± 3.77) and 
higher growth. In such a situation the only way in which the model 
can align itself with the measured fag is by choosing a smaller 
coupling parameter (^ 1 ). We have placed very stringent bounds on 
/(B) models with a ACDM bacgkround, constraining their only 
free parameter to be Bo < 1.36 x 10“^ (la C.L.). This is com- 


petitive with the constraint from solar system tests and clusters 

dHu & Sawickill2007l: 

ISchmidt et al. 20091) and other cosmologi- 

cal measurements ( Xul 

2015l;lDossett et al.l20l4lRaveri et al.l20l4 


IPlanck Collaboration et '^l2015b|). 


Finally, we have analyzed the growth index parametrization of 
the growth rate, measuring 7 = 0.612 ± 0.072 (11.7% measure¬ 
ment), which is completely consistent with the general relativity 
pr ediction. This is a sligh t improvement on the 16% measurement 
of ISamushia et al.l (120141) . We also note that our measurement of 
growth index is slightly less precise than current best measurem ent 
7 = 0.665 ± 0.0669 (10% measurement) I Johnson et al.l(l2015h us¬ 
ing combination of galaxy power spectrum, velocity power spec¬ 
trum, Type la SNe, the cosmic microwave background (CMB), 
CMB lensing, and the temperature-galaxy cross correlation. 

It is remarkable to notice that even after allowing many differ¬ 
ent kind of degree of freedom, our analysis shows that everything 
is consistent with vanilla ACDM cosmology and General theory of 
relativity. 


(c) 2015 RAS, MNRAS 000,[THT3] 
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